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Abstract 
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constraints on the correlation function from the transformation properties of the fundamental fields 
of QCD occurring in its definition. We further develop a method to construct an ansatz for this 
correlator. We present the complete leading order set of generalized parton distributions in terms 
of the amplitudes of the ansatz. Finally we conclude that the number of independent generalized 
parton helicity changing distributions is four. 
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I. INTRODUCTION 



In quantum field theory the non-perturbative nature of composite particles is described 
by matrix elements of operators of the fundamental fields of the theory evaluated between 
the initial and the final state of the particle under consideration. Accordingly, in QCD, the 
information on the internal partonic structure of hadrons is entirely given by matrix elements 
of all possible quark and gluon operators evaluated between hadronic states, forward and 
non-forward. 

Forward matrix elements of quark and gluon operators, i.e. operators between hadronic 
states of equal momenta, are investigated in hard inclusive processes, while non-forward 
matrix elements, i.e. operators between hadronic states of different momenta, have been 
investigated in recent years in the context of non-forward high-energy exclusive processes 
such as Compton scattering in the deeply virtual kinematical limit ( DVCS ) and hard 
diffractive vector-meson production. Experimental results Q, 0, 3, 0] indicate the feasibility 
of measuring these processes much more precisely with dedicated experiments in the future. 

In this article we will focus on hadronic matrix elements of bilinear quark operators, since 
these are the ones involved in the description of dominant contributions to non-forward hard 
processes. We consider their Fourier transforms, i.e. the quark-quark correlation functions, 
which contain all information concerning the non-perturbative nature of hadrons. From 
non-forward quark-quark correlation functions, for convenience, the so called generalized 



parton distributions (GPDs) are defined. 

Following and generalizing the method developed for the ordinary forward parton dis- 
tributions 3, 1^1 , we will formulate the most general ansatz for the quark-quark correlation 
functions starting from general principles, and we will then analyze leading order GPDs by 
projecting the ansatz with different Dirac matrices. In this way we will be able to express 
the GPDs in terms of the amplitudes entering the ansatz and establish a formal method for 
the determination and classification of the independent quark GPDs. 

The number of independent GPDs occurring in a given Dirac projection is not self- 
evident. For instance it has been debated in the literature 0, Q whether there are two or 
four independent quark helicity changing GPDs corresponding to the forward transversity 
distribution. With the method developed in this article we will show unambiguously that 
there are indeed four independent helicity changing GPDs at leading twist. 
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The outline of the work is as follows: in Section |H] we define the non-forward quark- 
quark correlation function and we introduce a new method to build an ansatz for it. We 
also derive constraints on the non-forward quark-quark correlator imposed by the hermiticity 
properties of the quark fields and their well-known behavior under parity and time reversal 
operations. In Section lllll we relate the non-forward correlation function to leading twist 
quark GPDs by tracing the ansatz with various Dirac matrices and by integrating over 
quark momentum components. Different Dirac structures probe different spin properties of 
the hadrons. Finally, in Section IIVI we draw conclusions and discuss the outlook for this 
subject. 



II. THE NON-FORWARD QUARK-QUARK CORRELATION FUNCTION 

We define the non-forward quark-quark correlation function $jA', j\(k, k', P, P'), depend- 
ing on the hadron and quark momenta (see Fig. for notation), by Fourier transforming the 
hadronic matrix elements of quark-quark operators 

^, i A(*,V l P,P / )= / ^^( fe+fc 'W 2 (P',A'|^(-i)^(|)iP,A). (1) 

We assume the hadron is a spin | particle, say a nucleonic target, which is in an eigenstate 
of light-cone helicity characterized by the initial and final light-cone helicity, A and A' 
respectively, which are defined from the spin vectors and S' 1 * 

S » = -(p»-™?- v '») s'^-iP'^-^-rv'n, (2) 
m K P+ ' ' m K P'+ ' ' V ' 

where m is the hadron mass, and v'^ is the null vector on the light-cone t>' M = [0, 1,0t]- 
Throughout the paper we use the component notation u = [u + ,u~,ut] with u^ 1 = (u° ± 
w 3 )/\/2 and the transverse part «y = (m 1 ,^ 2 ). Explicit representations of spinors for light- 
cone he 
Lepage 

element of quark-quark operators and not an expectation value like the forward quark-quark 
correlation function Q|. Therefore it provides more general information about hadrons 
compared to the ordinary quark-quark correlation function. 

The known properties of the quark-quark correlation function determine the structure of 
an ansatz, which will be the starting point for our analysis. The quark-quark correlation 



icity eigenstates are given for instance by Kogut and Soper [3| or Brodsky and 
1^ | . The non- forward correlation function in Eq. represents a transition matrix 
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FIG. 1: Diagrammatic representation of the correlation function $jA', jhik^k' , P, P'). From an 
initial hadron with momentum P a quark with momentum k is taken out, and reinserted with 
momentum k' = k + A to form the final state hadron with changed momentum P' = P + A. 

function in the helicity basis, denned in Eq. (0), is a 4 x 4 matrix in the partonic Dirac space, 
labelled by % and j, and a 2 x 2 matrix in the hadronic helicity space, labelled by A and A'. 
Therefore we formulate an ansatz given by the product of a partonic and a hadronic sector 
separately. The partonic sector is spanned by the set Tij of the 16 independent 4x4 partonic 
Dirac matrices, and the hadronic sector is represented by all possible independent spinorial 
products, Uk(P', A') T k i ui(P, A), evaluated between final and initial light-cone helicity states. 
Since the correlation function in Eq. (^) is a scalar in Lorentz space, we saturate the open 
indices of the tensorial product of hadronic and partonic sectors with all possible independent 
tensors, t fJil ... fJipU1 ... Uh , constructed from the kinematical variables k = (k + k')/2, P = (P + 
P')/2, and A = P' — P, the metric tensor g a/3 and one antisymmetric tensor e Q/ g po -, leading 
to the ansatz of the form 

$ W)jA {k, P, A) = f J 1 ' [u k (P', A') Y k r ^ Ul (P, A)] V ..,,,,,. ,,,(/>. k, A) . (3) 

In the hadronic sector the spinorial products (spinor indices suppressed) 

u(P',A') u(P, A) , u{P',A') 75 u(P, A) , u(P', A') a af3 u(P,A) , (4) 



u(P\ A') 7 Q u(P, A) , u(P', A') 7 °75 u(P, A) (5) 

occur. From the Gordon identities we can deduce that only the three spinor products in (J3J 
are independent. These spinorial products are evaluated in a specific frame in Appendix 1X1 
Furthermore, not all possible contractions of indices have to be taken into account in the 
construction of the ansatz for the correlation function. The relation [10] 

= u(P',A') u(P,A) A a + u(P',A') icr al3 2Ppu{P,A) (6) 
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shows that contractions of the tensorial spinor product u a a/3 u with Pp reduce to the scalar 
spinor product multiplied with A a , and the relation [l(3] 

= u(P', A') 75 u(P, A) 2P a + u(P', A') i a af3 7s A p u(P, A) (7) 

together with the known identities 

& at3 e a p p(7 = -2 ia pa 75 

v Q(3 e apaT = -i-f 5 {a pa gP - a pr g% + a aT g?) (8) 

entails 

u(P',A') a afi u(P, A) e a/3pa A p = Au{P',A') 75 u(P, A) P CT , 

«(P',AV"^(^A)e^ = u(P',A')Tfcu(P,A) {2P a gg. - 2P T g%) 

+ ±u{P',A') a^ s u(P,A) e l5aT ^ , (9) 

which reveals that the tensorial spinor product u a a/3 u contracted with either e al 3 pa A p 
or e apcrT A p is proportional to structures already accounted for. Note that the tensors 
t^— n P vi-v h {P 1 k, A) contain at most one Levi-Civita symbol since tensors with more than 
one do not result in new structures. The product of two Levi-Civita tensors reduces to the 
product of Kronecker symbols, from which no new structure arises. In a similar way the case 
of tensors t pi ... PpUl ... Uh (P, k, A) with more than two Levi-Civita symbols can be excluded. 
The correlation function fulfills the following constraints derived from properties of 



the Dirac quark fields and of the 
time reversal transformations f cf. 



radronic states under hermitian conjugation, parity and 



®l AljA (k,k',P,P') = 70$ iA,jA'(k', k, P', P) 70 (hermiticity constraint) 

^iA',jh(k, k', P, P') = 7o $i_ A 'j-A(fc, k', P, P') 7o (parity constraint) 

&*iA',jA(k, k', P, P') = (— ij^C) $jA'jA (k, k', P, P') (—ij^C) (time reversal constraint) , 



(10) 

where the notation u = (u°, —u) for momenta and spin vectors is used. These constraints 
can be implemented in building the ansatz for the correlation function. As far as the time 



reversal constraint is concerned, Collins has shown that the time reversal constraint is not 
applicable if a Wilson line is inserted in the quark correlation function. Therefore impo sing 
the time reversal constraint corresponds to defining the correlator without a Wilson line.|l2j] 
Taking into account all the above information leads to the most general ansatz for $. 
The lengthy expressions are an intermediate result of our investigation and are explicitly 
given in Appendix [Bj 



III. GENERALIZED PARTON DISTRIBUTIONS 

rpi 

The generalized parton distributions $^ A , jA (x, ^; t) are obtained as traces of the quark- 
quark correlation function with the Dirac matrices T, integrated over the transverse and 
minus components of the quark momentum, and k, respectively, 

$ l $ A (x,Z;t) = ~ J d 2 k T dk- Tr[$ A , A r], (11) 

where x and £ are the light-cone momentum fractions x = k + /P + and £ = — A + /(2P + ) and 
the Mandelstam variable t = A 2 denotes the momentum transfer squared. The projections 
$ A , A carry only hadron helicity indices since the parton Dirac indices have been saturated 
in taking the trace. Following j?[ 0| we rewrite the integral in Eq. (fTTjl with respect to the 
covariant integration variables, a = 2P ■ k and r = k 2 , and the azimuthal angle <fi 

4V^;*) = f dadrd<pe(xa-x 2 m 2 + ^-r)^0l. (12) 

The projections of the quark-quark correlation function $ with the various Dirac matrices T 
in Eq. (|T2*|) correspond to the different generalized distribution functions and determine also 
which distribution functions occur in a non- forward hard process at different orders in 1/ P + , 



P + scaling with the hard sea 



distribution functions 



BflQ 



e in the process. In particular the leading order ( twist 2 ) 
are obtained by projecting out the ansatz for the correlation 



function with the Dirac matrices 7 + , 7 + 75 and ia +l ^. 



A. Unpolarized parton distribution 

Let us consider the unpolarized generalized distribution functions of the proton. Thus 
we project the non-forward quark-quark operator with the matrix 7 + . Following ^| using 
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the notation of [18] the GPDs H(x,£;t) and E(x,£;t) are defined by 



*£a = \ f %-'" P "" <-P',A'l*(-^/2)7 + *( Z /2)|P,A)| I+ .„,r T =o 



u(P',AQ ? MP.A) t) + !(WW) B( () 

2P + 4m P + 



(13) 

Evaluating the quantity &2'A ^ or both proton helicity flip (A' = —A) and helicity non- 
flip (A' = A), the generalized distribution functions, H and E, result from the following set 
of two equations 







,h+] = _ 
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= V^T 2 (H-^E), 
$LV = - ( $!^1 J ) = 77 P (14) 



with t defined as 

4£ 2 7Tl 2 

and a phase 77 given by 



(15) 



A 1 + zA 2 



Substituting the expression of the ansatz (|B7|) and (jB8|) for the helicity non-flip and helicity 

flip correlation functions, $ ++ or $ , and $_ + or respectively, in (J12|) and tracing 

with 7 + we have a set of two equations for the two unknown distribution functions H and 
E, which is solvable and gives the two unpolarized generalized parton distribution functions 

as 



E = --^=A^(x^t) 



(17) 



where we have introduced the function A^(x,£;t) and A^ 3 \x, £; t) defined through the 
coefficients of Eqs. flB3) and (|B8l) 

A<U(x,€;t) = 1 / dadrd(j>6(x(T-x 2 m 2 + — - t) [djp + x Sp + 4£ 2 d$] (18) 



and 



x 2 t 



I (-20 [4 3> + *4 3) + <4? + 4^4?] 

_ 2? ^^_ 2 ,_e= 4 ,}. (19) 

The expressions for the generalized distribution functions H and E in Eq. (JT7j) do not contain 
any term in the amplitudes and dm , since these amplitudes are related to the amplitudes 
d$ and d { £\ respectively, as shown in Eqs. (jrl4j) . (fB5jl . (jrlfij) . Ml . (jBTTij) and fRTTj) . 

In the forward case rotational invariance around the collinear axis implies the conservation 
of the longitudinal component of angular momentum, i.e. it requires total helicity to be 
conserved (refer to FigJTJ) 

A + A' = A' + A . (20) 

Helicity conservation in Eq. (|20j) then shows a link between the quark and nucleon helic- 
ity degrees of freedom. Through the projection of the quark-quark correlation function 
with the matrix 7 + the nucleon helicity does not flip and the only possible hadron helicity 
combinations in the forward limit are + + and . 



B. Polarized parton distributions 

The polarized quark distributions, H(x,£;t) and E(x,£;t), are defined by the Fourier 
transforms of the axial vector matrix element 

*t + A 7sl = \ J %<?* P *'~~ <-P'.A'l*,(-^)7 + T 5 #'„(|)|i',A)U + =„,*.o 

= fi (P',A-hW.A) ~ H(Xt £; t) + .(P'.A^yP.A) S(x (; f) (21) 
For the different proton helicity combinations we now find 



[7 + 75] 



= (H_J? E) ^ (22) 



and 



$L^ ] = (^y) * = V £ E . (23) 
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Substituting the ansatze (jB7|) and (jB8|) for the correlation functions we obtain a set of two 
equations in the two unknown functions H and E 

H = - 7 =L= (BV>(x,t;t)- r J^ i i g (3 W;t) 

E = -— ^=^ 3 W;t), (24) 

where we have introduced the function B^(x,$,;t) and B^^x, defined as 



B w {x^;t) = 1 / dadrdcf) 9(xa-x 2 m 2 + — - r)[d$ + xd$ + ^ 2 d 
V' 1 < 2 ■' 1 



(25) 



and 



[di? + 4 e 2 4 3 8 } + x {d { s + 4 e 2 4?) + 4 e 2 <* 



2H 



+ ^= [2 « 4? - x 4?] - ^ [4fs + - 4?J r (26) 



respectively. 



C. Parton helicity changing distributions 

There are also twist 2 generalized distributions that change the helicity of the active par- 
ton. The corresponding quark distributions are constructed from the operator ip q i 0" + *75 i/j q . 

n 

By counting the helicity amplitudes Hoodbhoy and Ji [9] introduce two independent quark 
helicity changing generalized distributions corresponding to hadron helicity flip and non-flip 
terms. However, Diehl [10] claims that there are four independent parton helicity changing 
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generalized distributions denned by 



„■ +ia/3 A p _ 

u(P', A') ^ a +i 75 A) #r(z, £; *) + «(P', A') ^ A) # r (z, £; *) 

rrr 



+ u(P', A') - a7/3 «(P, A) E T (x, £; t) + u{P', A') & u(P, A) E T (x, £; t) 

2m m 



(27) 

In Eq. f)27|) the Lorentz index i takes the values 1,2). In the following we show that the 
same conclusion concerning the number of independent helicity changing GPDs is obtained 
by tracing the ansatz for the quark-quark correlation given in Eq. ()B7jl with the matrix 
ipqi <T + *75. Following Diehl we define 

^a = ^- (28) 

By comparing the definitions for Q l \>\ in Eq. (|28|1 and the projections of the non-forward 
correlation function $ with matrix i<j t+r )§ ) the GPDs, Ht, Ht, Et, and Et, arise. In fact 
the trace of the ansatz for the non-forward quark-quark correlation function with the matrix 
i < jl+ l5 gives a system of four linear independent equations 

$ltf 75] (*,£;*) = c^ + d^ 

$^ l+75l (x,e;t) = C w - 

$t i+75l (x,e;0 = C^* - Z^* (29) 
where C 1 - 1 ^, expressed in terms of the coefficients of the ansatz, reads 



I / T 2 t f p+VO" _ f pcr+i 

C {1)1 = ; / c/arfrc/0^(x ( x-x 2 m 2 + =— - r) z < ^-RA (7 4 1 ) + — =T^A CT d 

a/1 — £ 2 7 4 I roP+ mP + 



(30) 
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and reads 

= 1 J dadTd(j)6(xa-x 2 m 2 + ^-T) Ia i (£ a - 2 m 2 £ x + k ■ A) d 



+ l 1 ^ xa -2iT-xk.A)d§ + 2 -^{d { £-xd { £) }. (31) 

where k ■ A = x (m 2 — |) (£ + 1) — £ a — A; T • A T . 
The function C^ 3 -** can be written as 

x 2 1 I m 



if x t 
C (3) = / dadrdcj) 9(xa - x 2 m 2 + — r) 77 



n/^oV 1 -^ 



A ! A j 

+ — ? 4s (2 m 2 £ x - 2£k ■ P - k ■ A) -\ ? 2 £ 4o ( 2 £ A; • -P — 2£ P-xk-A) 

m 6 m A 



r-po+i pa+i A i 

-2 — *+ ^a, i 4? - 2 — g+ EpA, e 4? + — ] 



(32) 

and in terms of the amplitudes dm is 

L>( 3 )* = f dadrdcj) 9(xa - x 2 m 2 + — -r)ri \ , m ^ =[—2 % £ (rig + x df) ) 



+ — 2^d < i f j{x 2 m 2 -2xk- P + k 2 )} 



(33) 

Inserting in Eq. (J27)l the results for the spinorial products, given in the Appendix [XJ we 
obtain the four possible helicity combinations 

e+^PpA^^/T^e 



(l+£)m 



iH T 2me +3pr7 P p A a a/1 - £ 2 
2 (p+ A i £ + i e +^p p A (7 ) 



T ■ 



^- 1 P 1 A^ £ ^P £ K 



11 



(l + £)m 



H 7 

— i Ht 2m e +jfXT P p A a y/l - £ 2 
~ i P + A J + C, e + ^ pcr PpA a 

+ / 



(35) 



QL 



2H 



Hn 



-iP+Ai + e + jP a P p A a 
e+jprpA (4m 2 £ 2 -^t + t) 



+ 4mE T £, 



, -P + A j + i e +jpCT P p A a 



-2H 



_P+ a-?' + i e + ^P p A ff ) 
iP+A^ + e +jpa P p A a 



Hn 



+ 4m£ 2 £ r 
- 2miEr 



P+Ai + i e +jpa P A c 



P + A j + i e+i^PpAe 



(36) 



(37) 



Since we are only interested in the number of independent GPDs we refrain from isolating 
H T , H t , E t , and E T in Eqs. (JMD, (J3SJ), © and (J37|). Instead, from the fact that Eq. 1(25) 
constitutes a set of four linearly independent equations, we conclude that one can define four 
independent GPDs from it. Thus, we have shown that the number of independent helicity 
changing generalized parton distributions is four, as claimed by Diehl [lflj ]. 



IV. CONCLUSIONS 

We have presented a detailed analysis of the non-forward quark-quark correlation func- 
tion. Constraints on the correlation function, not yet in the literature, were obtained by 



12 



implementing the known properties of the fundamental fields of QCD, quarks and gluons, 
under parity and time reversal transformations and applying hermiticity. We developed a 
new method to construct an ansatz for the correlation function. The quark-quark correlation 
function could then be expressed in terms of tensorial structures formed by the independent 
dynamical vectors and by Dirac matrices. The constraints obtained were implemented to 
reduce the number of independent amplitudes multiplying these tensorial structures in the 
ansatz. 

Finally we projected out the leading order GPDs, i.e. we expressed the unpolarized, 
polarized and parton helicity flip distributions in terms of the amplitudes occurring in the 
ansatz. The formalism adopted allowed us to conclude that the number of independent 
parton helicity changing distributions is four, in agreement with Diehl's argument [10|. We 
stress that the result about the number of the independent GPDs was obtained by Diehl 
in a completely different way and this is a confirmation of both methods used to approach 
the problem. On one hand we wrote the most general ansatz which can describe non- 
forward quark-quark correlation functions, i.e. we represented matrix elements of non-local 
non-forward quark-quark operators in terms of tensorial structures, built from the involved 
momenta on the basis of general properties of invariance. Then we traced the ansatz for the 
non-forward correlation function with different Dirac matrices and we could read off which 
of these structures contribute to each GPD. On the other hand Diehl's approach was to 
count the number of independent helicity amplitudes occurring in DVCS cross sections on 
the basis of time reversal and parity invariance which these amplitudes have to fulfill. 

The advantage of having built an ansatz for the non-forward quark-quark correlation 
function is that, by tracing it with the different T Dirac matrices, we gain the different 
generalized distribution functions in terms of some of the amplitudes occurring in the ansatz, 
and we are thus able to predict the dependence of GPDs upon the different fundamental 
structures entering the ansatz. 

Note that in Eq. (fTTJ) we integrate over d 2 d k~ and thus we consider only distribution 
functions which do not depend on the transverse momentum of the quarks kx- We remark 
that, in principle by having an ansatz for the off-forward quark-quark correlator, one could 
extract generalized profile functions which depend additionally on the transverse momentum 
kr of quarks. For instance the investigation of ky-depending ordinary parton distributions 
has been extensively carried out by many groups theoretically and experimental investiga- 
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tions are currently under way. For non-forward processes no formalism for the systematic 
study of kr-dependence has yet been attempted and an experimental program on effects 
seems far beyond present abilities. In the foreseeable future there are good prospects to 
acquire some knowledge on GPDs depending on (x, £, t; Q 2 ), as well as possibly additional 
k T - dependence. 

We have worked out a powerful method of analysis which in the present paper was applied 
completely to the leading twist level. The same method can be implemented to investigate 
twist 3 and twist 4 generalized distribution functions. For instance one could expect that 
useful relations between leading and next to leading order generalized distributions could 
emerge as suggested by similar experience in the forward case. In this sense the present 
work represents a valuable starting point for further investigations. 



APPENDIX A: SPINORIAL PRODUCTS 



The independent spinorial products of ansatz J3J) listed in Ec 



by using explicit expressions for light-cone helicity spinors 

/ V2P+ + m\ 
1 P 1 + iP 2 

(2 v / 2P+) 1 /2 V2P+-m 

\ P 1 + iP 2 ) 



14] 



U LC (P, +) 



u L c(P, 



are most easily evaluated 



/ -P^iP 2 \ 

V2P+ + m 

P l -iP 2 

\ -V2P + + m ) 
(Al) 



{2V2P+y/ 2 



normalized according to 



u{P, +) u{P, +) = u{P, -) u{P, -) = 2 m 

u{P, -) u(P, +) = u(P, +) u(P, -) = . (A2) 

In a frame of reference where the longitudinal direction is defined by the proton average 
momentum P, the momenta of incoming and outgoing hadrons are parameterized as (in 



the light-cone component notation 
two-dimensional transverse vector) 

P' 1 = (P-A/2f : 

p'f = (P + A/2f : 



\z + , z , z±], where z d 



± z 3 )/y/2 and z±_ is 



m 2 + A^/4 Aj 



2(1 + 0P + ' 2 

m 2 + A^/4 A ± 
2(1 -0P + ' +_ 2" 



(A3) 
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With this parametrization one obtains for the spinorial products in the helicity non-flip case, 
i.e. A' = A = ±1 



u{P',A') u(P,A) 



2m 



u(P',A') l5 u(P,A) = A 



2m£ 



1 -£ 2 



u(P',A') a^ v u{P,A) 



i m (AV" - A v v' p ) t A 2m e pupa P p v' a 



+ A — 



The spinorial products if the helicity is flipped and A' = —A = 1 are 



u(P',A') u(P,A) 
u(P',A') l5 u(P,A) 
u(P',A') a^u(P,A) 



-ri Ay/t -t 

2 i (P P A U - P u A tl ) 

y/t^t 



2 i m 2 (AV - A u v' p ) 



P+(l-e)Vh=t 



(A4) 



while if A' 



-A 



4Am 2 £e pupa P p v' a 2A t pupa P p A a 



P+{l-e)Vh=t y/t^t 



(A5) 



u(P',A') u(P,A) 
u(P',A') l5 u(P,A) 
u(P',M) u{P,A) 



V*AVt^t 

2i (P P A U - P U A P ) 



2 % m 2 (A p v' u - A V) 
P+(l-e)y/h=t 



4Am 2 £e pupa P p v' a 2A e pupa P p A a 



P+(l-e)Vh=t y/t^t 



where the phase factor is given as 



Ai + iA 2 
IA, I 



and 



IA 



-4£ 2 m 2 4£ 2 m 2 + A 2 



+ 



(A6) 



(A7) 



(A8) 



i-e i-e 

which contains the implicit definition of the quantity to- Studying the form factor decom- 
position of the tensor current of the proton [lfj, we need additionally the following Dirac 
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bilinears. For the helicity non-flip case, i. e., A' = A = ±1, we have 

-i P+ AA j + £ e +jpa P p A a 



u(P',A') e+^P pla u(P,A) 



^2 



U 



(P',A') a +j 7 5 u(P, A) 



u{P', A') e +jpa P p A a u{P, A) = -2mie +jpa P p A a y/l - £ 2 



The results for the helicity flip case with A' = — A = +1 are 



u(P',A') e +J " ff P p7ff tt(P,A) = -77 

u(P',A') a +j 7 5 u(P, A) = -77 

u(P',A') e +jpa P p A a u(P,A) = -r] 

u(P',A')e +jpa A pla u(P,A) = - v 
and for A' = -A = -1 

u(P',A')e^^P pla u(P,A) = if 

u(P',A') ct+V u(P,A) = T]* 
u(P',A') e +jpa P p A a u(P,A) = rf 



m 



2 £(P+ A^'A + i e +jpa P p A a ) 



4m£ 



iP + AA j + e +J > a P p A CT - 
e +jpa P p A a (4m 2 £ 2 -£ 2 t + t) 

2 P+AA- 7 ' + i e+^PpA,, 



- 2 



m 2 £(P+ A* A + i e +jpu P p A a ) 
2 iP + AAJ +e +jpa P p A a - 
e +ip°p p A a (4m 2 £ 2 - £ 2 1 + t) 



a 



(P',A')e+^A p7(jM (P,A) = rf 



4m£' 



P+AAi + i e +jpa P A 



(A9) 



(A10) 



(All) 



where the phase is again given as in Eq. (jA7|) . Note that latin indices i,j = 1, 2, while greek 
indices run from to 3 and e +1 ~ 2 = 1, since we have adopted the convention e 0123 = 1. 



APPENDIX B: ANSATZ FOR THE NON-FORWARD QUARK-QUARK COR- 
RELATION FUNCTION 

For the sake of convenience, in view of the numerous coefficients, we now switch over from 
explicitly displaying the hadron helicity indices to a more dense notation using the index 
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(k) defined as (1) = ++, (2) = H — , (3) = — h, and (4) = and suppressing spinor 

indices ij. Moreover contractions of momenta with Levi-Civita tensors will be written in 
the compact notation e al3lS ciabpc^ds = e abcd . 

In order to write the ansatz for the quark-quark correlation function as indicated in 
Eq. (J3J) we now multiply the independent spinorial products, given in Eqs. (jA4j) . (jA5|) and 
(IA6|) . with the 16 independent partonic Dirac matrices and saturate the free indices with 
the following tensors ( indices a and (3 are intended saturated with spinorial products, while 
indices /i, v with indices from partonic Dirac matrices ): 

1 P^ e ^PkA P^k v P^A U k^A U e ^Pk e fj,uPA e /ii/kA 

g^k 13 g afJi A 13 A a WP» A a Wk» A a WA» ^ af3p e^ af3 ~ k 

e aj3Pk pfi e af3Pk ^ ^afiPk ^ji ^fiPk ^fiPk £J3 ^PkA £J3 g ^ ^fiPkA 

9/xa 9u/3 9{ia 

¥P V g m ¥k v g m ¥A v g m ^P v g m A p k v g m A 13 A v 
k a A 13 P^k u k a A 13 P^A V PA^FA" e^ af3 

e a/3Pk pp, J,u e a/3Pk p/j, ^u ^afiPk j^/i ^u 

^fia/3P pu ^fiaPPj^u ^a(3P ^u ^ia0k pu ^afSkJ^u ^ictf)k 

{V-voiP^fi e fj,vaP £J3 e fiuak TS ^uak ^vPk j^a ^(3 

g m e^ g, a e» pkA ¥ g m ^A? 

^iivPA k a A^ e /J-vkA k a A@ ^Pk fofl pu e afiPk Jj/ 

e a^Pk fo(3 ^u ^afiPk ^ pu e a/iPk ^(3 jj> £ afiPk £J3 j\u 

e !iPkA yx £J3 pu e MPfcA£a A /3^ ^PkA j^a ^u ^ 

The method introduced in Eq. (JHJ) produces an ansatz containing Lorentz tensorial structures 
multiplied with various amplitudes that are functions of all possible Lorentz scalars. In the 
ansatz we choose to indicate explicitly only the dependence on scalar products A-v ' involving 
the vector v' which occurs in the definition of the helicity eigenstates of the hadrons. Through 
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a renaming of the amplitudes each tensorial structure is thus multiplied by an amplitude 

(k) 

dm depending implicitly on all scalar products which can be formed from the momentum 
vectors k, P and A 



M). 

For the helicity non-flipped case with k — 1,4 the ansatz reads 
<&<">( jfe.P, A) = 



(B2) 



+ 



■ A-v' 



(K ) mA-v' ( K) 
ma\ + " 



_4 K) + ^4 K) + ^4 K) + ^ ) + ^ A ' w, 



p + "10 +Fp + "11 



a\ 7 > +ri — a\ 8 > + a\ 9 > + -£= =— 

'A" ^ a iat P*A v A • r „, a^A^fc" ™ <J IIV A^ 
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+ 7 5 faff + 7 5 f ^ a<? + 7 5 ? a£> + 7 5 ? ^ 4? + 7 5 Aa W 
5 

+ 7 5 A' 



(«) 



p+ "20 

a^P^A" (K) <vP"A" A • u' (K) a^A^fc" (k) <vA"fc" A • 

+ °21 + fTT - a 22 + fl 23 + ~ 

21 m P+ m - 1 



p+ tt24 



p+ °26 + cr /i//^ ^ p + fl 27 + m << 



'28 



1> 29 + 
-fa^AW (K) , -fa^AW A-v' 

+ «32 + " 

m m P + 



7 5 yP"fc" A 
m 

„5 



t^p^a; a • w 

'30 + „„ p + a 31 



'32 

m 



^ p+ fl 35 + 7 V n U p+ "36 + m 2 



a (, ° + V 5 a P»v' v —a { ' 

p + "33 + 7 CT/i^-r V p + a 3 

; ,„ (*) , 7m ^ A A- 

p + "36 + """ - ' ; '~ " ' ' - 



V' 



m 

7mJ_ _(«) , 7m ^ («) 7m g (ic) 7 7m ^ («) 

p> "39 + p> "40 + p> "41 + ^2 



p+ fl 38 



+ >y 5 a^v* 

>\ "I" ^ "ill ^ 5 u 42 

7 5 7m^" A A • («) ^ 7 5 7m^'^ (ic) , 7 5 7m^ PA ^ W , 7 5 7^^ A "' (k 

+ m 2 p+ fl 43 + p+ "44 + p + "45 + p + °46 



,(«) 



+ 



e PkAv' 



) ( K ) i 
7TT P+ 



^4? + H$ + + ^4? + 7 5 f 4? 



"55 + 



"56 "I 



m 



l 57 



(B3) 



For the helicity non-flipped case we implement the constraints on the non-forward quark- 



18 



quark correlator imposed by the hermiticity properties of the quark fields and their well- 
known behavior under parity and time reversal operations as stated in Eq. (|10p. In particular 
parity invariance imposes the following relations 

a K m = a K m m = 1, 2, 5, 6, 7, 8, 9, 10, 11, 19, 20, 21, 22, 23, 

24, 25, 26, 27, 42, 43, 44, 45, 46, 48, 52, 53, 54 
a K m = -a* m = 3, 4, 12, 13, 14, 15, 16, 17, 18, 28, 29, 30, 31, 32, 33, 34 
35, 36, 37, 38, 39, 40, 41, 47, 49, 50, 51, 55, 56, 57 



(B4) 

Assuming hermiticity 

a K m = (a*)* m = 1, 4, 5, 6, 7, 10, 11, 12, 14, 17, 18, 19, 22, 24, 25, 26, 

29, 30, 32, 36, 37, 39, 42, 45, 46, 48, 51, 54, 56, 57 
a n m = -«J* m = 2,3, 6, 8, 9, 13, 15, 16, 20, 21, 22, 23, 27, 28, 31, 

33, 34, 35, 38, 40, 41, 43, 44, 47, 49, 50, 52, 53, 55 



(B5) 

Imposing the time reversal constraint reduces the number of independent amplitudes in the 



ansatz since 



< = «J* m = 1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 

13,14,15,16,17,18,55,56,57 
a K m = -«J* m = 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 

31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42 
43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54 



(B6) 

From Eqs. (jB4|) . (jB5|) and (jB6|) the diagonal amplitudes are either real or pure imaginary. 
For the case in which the hadron helicity is conserved, ( k = 1,4 ), the ansatz for the 
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quark-quark correlation function can be defined as real by multiplying the pure imaginary 
amplitudes by the imaginary unity i 



k,P,A) 



( K ) i 5 

m a\ ' + 7 



a\ + fa^ + f a? + A 



P+ 



10 



» 



A-v' (<s 

,2 



p + tt 17 + 7 f p + "18 



+ i 4? + i 4? + i 4? + i vAve 4?) 

m P + m m P + 



m 



a 28 ^ 



m 



P+ 



a 31 + i 



m 



P+ 



ag 7 + z 



P+ 



' 10 



+ i 



7m ^ W , 7 5 7m^ A W , 7 5 7^ PzW W , 7 5 7m^' („) 



P+ 



a 41 + 



rrr 



«42 - 



p+ 



a 45 + 



PkAv _ _ 

- i ^rr I m a$ + fag + % a { £ + 7 5 f 4? + 7 5 ? a$ + ^ 



<r„„P^ u 



m 2 P^ 



,(«) 

'56 



: a 57 
z m 



(B7) 

If the hadron helicity is flipped from Eqs. (|A5|) and (|A6|) one deduces that the ansatz has 
a factor — rj (77*), given in Eq. ()16j) . which reflects the difference in phase of the initial and 
final hadronic spin. From the method introduced in Eq. (J3J) the ansatz for the off-diagonal 
components of the quark-quark correlation function ( k = 2, 3 ) reads 



*W(jfc,P,A) - ^-^3) 



(«) mA V (K) 5 (K) 5 mA V (k) (k) 

ma l H p> a 2 +"^7 a 3 +7 pTf a 4 + f a 5 

2 a ../ a 2 



+ + fi^pf) 4" + P«. M + + 4? + £4? 

+ ^4? + ^)4f + 4? + ^ 4? + 7 5 f 4? 
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• ff^pf 4? + 7 5 f (^) 2 4? + 7 5 P 4? + 7* ^ 4? + 7 5 P 4? 

+ 7 5 44? + 7 5 4^ 4? + 7 5 U^pf) 4? + 7 5 4? 

5 AVm 2 w <yP>fc" (K) a^P^A-v' {k) a^Mz" A ■ v' 2 (k) 
P+ P+ m m P+ m P+ 

<Vp>A" {K) a^A" A - v> w yM- A • v> 2 (w) cr^ v k^A u (<s) 
+ 030 + ^ P^° 31 + ^ ° 32 + m 33 

^>A- A • ,/ (K) ^>A- A • „' 2 w m w 

+ ^pv^ ( ^) 2 4? + ^sv^ 4? + a^jk (^pf) 2 4? 
+ vav5 4? + ^ ^ 4? + 7 5 ^ 4? 

i? + 7 5 ^(^Mi? + 7 5 ^4f 
m P+ m 

3 V + 7 (-pr-) 4 7 J + 7 ~ ay 

# + 7 5 o, v PN^ og + 7 S a^PV^ (^) a<? 

+ 7 5 a^V^ 4? + 7 5 (^) 4? + 7 5 <vAV£ 4? 

5 w m A-v' {K) 7/ ,€^ A (w) y^f^ A-v' ((c) 7^^ (k) 

7 ,e^' A • „' (K) 7fi e^ (K) 7m A • */ {k) 7fx e^ w 

P+ P+ P+ P+ P+ P+ 

+ 7 m 2 a 63 +7 m2 p + a 64 +7 p+ a 65 +7 p + p + a 66 

4 ~5 7 ^ («) , 5 7m ^ (ic) 5 7m ^ («) , i r M 

A-^' («) . s («) , m^' v> («) , « («) , « («) , /v («) , /v A • v' (K ) 
+ -p^- a n + r a 72 + f-p^a 73 ; + fta y 74 > + ^a 75 ; + 4ka 76 ; + 4~p^ a 77 

_L ^5 J K ) 1 -,5 A ' V' («) 5 p («) 5 p AV ( K ) 5 -j ( K ) 

+ 7 a 78 + 7 p + a 79 + 7 f a 80 + 7 f ^ + a 81 + 7 f a 82 





A-v' 


m 


P+ 


o^A" 


A-v' 


m 


P+ 


o^A v 


A-v' 


m 


P+ 
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A V 



+ 7 5 ^ ^— a £> + 7 5 Aa^ + 7 5 A^ C + 



(«) 



A V 



(«) 



P+ 



P+ 



(k) OjiV 
u 86 



a„ v P»k u A ■ v' 



P+ 



.(«) 

l 87 



+ 



tyP^A" {K) a^A^A-v' (w) | (y^A" w , y^A-D' (k 



ass + 



to" 



P+ 



^89 + 



to" 



a 90 H 



to" 



P+ 



"91 



+ m , 



mag +m 7 5 4 K 3 ) + HI + + +7 5 f- (K) 

+ 7 5 Aag + ^^a^ 
to 



■a^ + 7 5 pa. 



98 



J-^ (K) tyi^A" (K) cy^A" ( K) 5 (K) 
to ° 100+ to fll01+ to 0102 + 7 to 0103 

5 ^A- w a^A 
-r 7 a 104 -t- 7 

to m 



c riPkA ^uPkA 

. >) + 2^1 >) + ^5 7^ («) 

a 105 + _ 2 "106 + 7 2 "107 



to a<£ + m 7 5 < 09 + f a& + # + MS + 7 5 f a& 



,5 » 



» 



(«) 



.,5 S («) , 5 /V («) , <V P ^ (*) , G p.vP iX ^ V («) , <yfc"A" (K) 



+ 7 5 ^iii + 7 5 « + 



+ 7 z «ii9 + 7 



to 
m 



aiie + 



to 



a\i7 + 



to 



l 118 



c fj,PkA 



cliPkA 



tt 120 + 7 tt 121 + „ 9 "122+7 - tt T5- 



to 



to" 



to" 



123 



(B8) 

We implement again hermiticity, parity and time reversal invariance. In particular parity 
invariance imposes the following relations 



d 2 m = V 2 d 



2 ^3 
rn 



to = 1,2,5,6,7,8,9,10,11,12,13,14,15,27,28,29,30,31, 
32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 63, 64, 65, 66, 67, 68, 
69, 78, 79, 80, 81, 82, 83, 84, 85, 92, 94, 95, 96, 100, 101, 
102, 107, 108, 110, 111, 112, 116, 117, 118, 123 

to = 3,4,16,17,18,19,20,21,22,22,23,24,25,26,42,43, 
44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 
60, 61, 62, 70, 71, 72, 73, 74, 75, 76, 77, 86, 87, 88, 89, 90, 
91, 93, 97, 98, 99, 103, 104, 105, 106, 109, 113, 114, 115, 119, 
120, 121, 122, . 
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(B9) 



Applying hermiticity 

d 2 m = (d 3 J* m = 1,4, 5, 7, 8, 10, 12, 14, 16, 18, 19, 21, 23, 25, 27, 29, 31, 

34, 36, 37, 38, 39, 41, 43, 45, 47, 48, 51, 53, 54, 57, 58, 61, 
64, 65, 68, 73, 75, 76, 78, 79, 81, 83, 84, 87, 88, 90, 92, 94, 
95, 97, 98, 100, 104, 105, 108, 110, 111, 113, 114, 116, 120, 121 

d 2 m = -(d 3 J* m = 2,3, 6, 9, 11, 13, 15, 17, 20, 22, 24, 26, 28, 30, 32, 33, 

35, 40, 42, 44, 46, 49, 50, 52, 55, 56, 59, 60, 62, 63, 66, 67, 
69, 70, 71, 72, 74, 77, 80, 82, 85, 86, 89, 91, 93, 96, 99, 101, 
102, 103, 106, 107, 109, 112, 115, 117, 118, 119, 122, 123 . 



(BIO) 

Imposing the time reversal constraint reduces the number of independent amplitudes in the 
ansatz since 

d K m = rj 2 m = 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 

42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 
58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 
75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 100, 101, 102, 103, 
104, 105, 106, 107, 116, 117, 118, 119, 120, 121, 122, 123 
d K m = -rj 2 (d 1 ^)* m = 1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 

19, 20, 21, 22, 23, 24, 25, 26, 86, 87, 88, 89, 90, 91, 92, 93, 94, 
95, 96, 97, 98, 99, 108, 109, 110, 111, 112, 113, 114, 115, 



(Bll) 

Because of the relations in Eq. (jB9|) . (|B10|) and (jBlljl some of the amplitudes are zero. We 
refrain from rewriting the lengthy expression of the off-diagonal ansatz. 
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In the forward limit A = the spinorial products in Eq. (jA5|) and Eq. (|A6|) cannot be 
built through the trace method since the product u(P,—)u(P,+) = u(P,+)u(P,—) = 0. 
This implies that the off-diagonal part of the ansazt does not converge. The method here 
developed is therefore applicable only to the cases for which A is different from zero. 
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